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Abstract 

Best-fit values of recent global analyzes of neutrino data imply large solar neutrino 
mixing, vanishing Ues and a non-maximal atmospheric neutrino mixing angle 023- 
We show that these values emerge naturally by the hypothesis of "scaling" in the 
Majorana neutrino mass matrix, which states that the ratios of its elements are 
equal. It also predicts an inverted hierarchy for the neutrino masses. We point out 
several advantages and distinguishing tests of the scaling hypothesis compared to the 
Lf, — Lfj^ — Lr flavor symmetry, which is usually assumed to provide an understanding 
of the inverted hierarchy. Scenarios which have initially vanishing Ues and maximal 
atmospheric neutrino mixing are shown to be unlikely to lead to non-maximal 623 
while keeping simultaneously Ues zero. We find a peculiar ratio of the branching ratios 
^ — > 67 and r ^ 67 in supersymmetric seesaw frameworks, which only depends 
on atmospheric neutrino mixing and results in r — > 67 being unobservable. The 
consequences of the scaling hypothesis for high energy astrophysical neutrinos at 
neutrino telescopes are also investigated. Then we analyze a seesaw model based on 
the discrete symmetry D4 x Z2 leading to scaling in the low energy mass matrix and 
being capable of generating the baryon asymmetry of the Universe via leptogenesis. 
The relevant CP phase is identical to the low energy Majorana phase and successful 
leptogenesis requires an effective mass for neutrinoless double beta decay larger than 
0.045 eV. 
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1 Introduction 



Observed lepton mixings are consequences of a non-trivial structure of the neutrino mass 
matrix Ai^- This symmetric matrix for Majorana neutrinos (having entries 771^3 with 
a, (3 = e, fi, t) is in the charged lepton basis diagonalized by the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) neutrino mixing matrix U. The very different structure of U compared to 
the quark sector for all possible neutrino mass orderings is indicative of an unexpected tex- 
ture of the mass matrix, and could hold important clues to our understanding of the physics 
of fundamental constituents of matter. To unravel this new physics, various Ansatze for 
Aiu have been made in the literature [1] and their associated symmetries have been sought 
after. One particular proposal, recently proposed by two of us (R.N.M. and W.R.), on 
which we will focus in this note, is called "scaling" |2j. The scaling hypothesis demands 
that the ratio is independent of the flavor a: 

"^T^ _ _ rrir£ _ ^ fixed j3 and 7 . (1) 

There are three possibilities and the only one phenomenologically allowed is when (3 = n 
and 7 = r. We shall call this case scaling henceforth. The resulting mass matrix reads 



(2) 



Similar matrices have been found independently in the context of specific models (see 
Ref. [3|). The most important phenomeno logical prediction of scaling is that Eq. ([2]) leads 
to an inverted hierarchy with rris = and Ues = 0. Atmospheric neutrino mixing is 
governed by the "scaling factor" c via tan^^23 = i-e., is in general non-maximal 

because c is naturally of order, but not equal to, one. It is interesting to note that current 
data analyzes (though at the present stage statistically not very significant) yield non- 
maximal tan^ 6'23 = 0.89 as the best-fit point [1] (see also [5], where the best-fit value is 
tan^ 623 = 0.82). The reason is that in the SuperKamiokande experiment there is an excess 
of sub-GeV electron events, but no excess either of sub-GeV muon events or of multi-GeV 
electrons. In a realistic 3- flavor analysis this prefers cos 6*23 > sin 6*23 [HE]. The best-fit 
value of Ue3 is zero. 

If these values, ^^23 7^ 7r/4 and Ues = 0, are indeed confirmed by future data then one 
should look for symmetries and/or models which are capable of predicting such a situation. 
Ideally, such a candidate should be rather insensitive to radiative corrections and should 
not require much, if any, breaking to achieve the values sought for. Scaling is one such 
appealing possibility, several general aspects of which will be discussed in Section [2l Typical 
models are shown to predict - when constructed in a supersymmetric seesaw framework 
- a characteristic ratio of the branching ratios fi —>■ ej and r — >■ ej. The latter decay is 
then too rare to be observable in presently foreseen experiments. Simple phenomenology of 
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fluxes of liigli energy astropliysical neutrinos at neutrino telescopes is predicted and studied 
in Section [31 We argue furtlier in Section H] tliat it is difficult to obtain 6^23 7^ vr/4 and 
Ue3 = in scenarios in which initially 623 = vr/4 and Ue3 = holds. Stressing that scaling 
predicts an inverted hierarchy leads us to compare the Ansatz with the flavor symmetry 
Le — Ln — Lt [6] . The latter is usually assumed to be the origin of an inverted hierarchy. 
We show in Section [5] that scaling possesses several advantages over L^ — L^ — L^-. If future 
experiments indeed show that neutrinos obey an inverted hierarchy, then one needs a full 
list of possible scenarios that can predict it. Necessarily, these are unusual symmetries as 
typical GUT models do not lead to an inverted ordering. Accordingly, we investigate in 
Section [6] a seesaw model leading to scaling based on the discrete symmetry D4 x Z2. We 
show that the baryon asymmetry of the Universe via leptogenesis can be reproduced and 
analyze the connection to the low energy parameters. We summarize in Section [71 



2 General Properties of Scaling 

The neutrino mass matrix giving rise to scaling is given in Eq. ([2|). It appears in the 
Lagrangian 

C = ^l^M,UL+'GMiiL , (3) 

where Ai^ is the charged lepton mass matrix. Diagonalizing it via Me = ViMf'^'^Uj and 
the neutrino mass matrix with U; Mf""^ Ul = M^, gives us the PMNS matrix 



U 



C12 Ci3 



"■S12 C23 — C12 S23 Si3 e 



\ ■5l2 523 ~" C12 C23 Si3 



-312 Ci3 



C12 C23 ~ S12 S23 Si3 e 

^Ci2 S23 — Si2 C23 Si3 e* 



si3 e-'^ \ 

S23 Cl3 
C23 Cl3 / 



p 



(4) 



where c 



cosOij, Sij = sinOij, and P = diag(l, e*"", e**-^^^-*) contains the Majorana phases 



[Tj. The best-fit values as well as the allowed la and 3cr ranges of the oscillation parameters 
are [Ij: 
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0.3, 0.8 r 10"' eV^, 



10"=^ eV 



(5) 



where Ami, = — mf and Am\ = \m'^ — ml\ 



We will most of the time assume that scaling holds in the charged lepton basis, i.e.. Up = 1. 
We nevertheless stress here the following interesting property of scaling: if Ui is non-trivial, 
then in the charged lepton basis we have ^Ay = Uj M.y Ue. In case Ue is only given by a 
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23-rotation it is easy to show that the ee entry of A^i, is not affected, and in addition we 
have: 

_ fh^^ _ _ c cos 6*23 - sin Oj^ ^ „ 

fner f^rr COS 622, + C Sm 622, 

where ^23 is the rotation angle of and m^^ are the entries of M.^. Consequently, the 
texture of the neutrino mass matrix is not changed at all, it still obeys scaling and predicts 
ma = ?7e3 = but now with tan^ 6'23 = 1/c^. If for some reason in M.y the scaling factor 
c is much larger or smaller than 1, then a 23-rotation from the charged lepton sector can 
still save the scaling hypothesis. For c ^ 1 we find tan^ ^23 — tan^ 6'23, while for c <^ 1 we 
have tan^6'23 — cot^6'23. This observation can simplify the construction of models. 

Interestingly, the characteristic predictions = Ue3 = are not subject to any radiative 
corrections when going from high scale down to low scale. This can be understood by 
letting RG effects directly modify the mass matrix, as done in Ref. [2], or by glancing at 
the RG equations of ^13 and m^. For both of them it holds that [S] 

772,3 oc m3 and ^13 oc 1713 , (7) 

i.e., 6'i3 = m3 = is stable under RG evolution. 

What are other phenomenological predictions of scaling? First of all, there will be no CP 
violation in oscillation experiments because of ^13 = 0. Then we note that from ms = 
it follows that one Majorana phase is unphysical. The other one appears in the effective 
mass to which neutrinoless double beta decay is sensitive [9J . This parameter takes a very 
simple form for the inverted hierarchy with = 613 = 0: 

(m) = v^Ami \/l - sin^ 2^12 sin^ a . (8) 

The range of (m) lies for the best-fit parameters from Eq. ([5]) between 0.019 and 0.051 eV, 
while at l(3)o" it ranges between 0.017 and 0.053 eV (0.011 and 0.057 eV). If the parameters 
do not conspire to render (m) at the very low end of this range then next-generation ex- 
periments will definitely observe neutrinoless double beta decay |9]. The conditions under 
which one can extract a from an observation of neutrinoless double beta decay are given in 
Ref. pro]. Unlike many other approaches, the scaling Ansatz can therefore be completely 
reconstructed. Of course, the scaling Ansatz is easier to disprove than to prove. However, 
if future experiments give very strong limits on {Uesl inverted hierarchy is present, 

then this would be a very strong hint towards the realization of scaling. 



How can a low energy mass matrix like Eq. ([2]) be achieved? We work of course in the 
framework of the seesaw mechanism ^JJ in which Ai^ = ^d, where is the 

Dirac and Mr the heavy Majorana neutrino mass matrix. One remarkable property of 
scaling is the following: if the Dirac mass matrix obeys scaling, i.e.. 



D 




(9) 
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then Aiu takes the form obeying scahng^ from Eq. ([2]) regardless of the structure of Mr\ 
Note that is not necessarily required to be symmetric. 

Within supersymmetrized seesaw models one has an interesting connection to lepton flavor 
violating (LFV) decays of charged leptons such as /i — > 67 [12]. RG evolution within 
theories of universal boundary (mSUGRA) conditions leads to off-diagonal entries in the 
slepton mass matrix, which trigger effects of LFV. In particular, branching ratios of the 
decays ii — > £^7 with = (/x, e), (r, e), (t, /i) are given by [12] 



a 



(3 + a, 



1)^1 



where (L)j 



Sij In 



{Ml)LMD)rj tan^/? 



Mx 
Mi 



(10) 



Here fwk = 174 GeV, Mj are the heavy Majorana neutrino masses and Mx > Mi is the 
scale at which universal boundary conditions are implemented. The SUSY parameters are 
Oq = Ao/mo with ttlq the universal scalar mass, Aq the universal trilinear coupling, and 
ms is a typical SUSY mass. Neglecting the only logarithmic dependence on the heavy 
Majorana masses, the branching ratios are proportional to the modulus-squared of the 
off-diagonal entries of M|, M^ . It follows with Eq. that 
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BR(r ^ 67) 



BR(r ^ez/z7) BR(/i ^ 67) 



{MIMd)31 



{MIMd)21 



^ = tan^ 6'23 , 



with BR(r evV) ~ 0.1784. The two branching ratios are therefore simply related by 
the atmospheric neutrino mixing angle. Consequently, such models predict them within 
(see Eq. (E])) a factor of less than four. The current limit of BR(/i cj) < 1.2 x 10~^^ [T3] . 
implies therefore that BR(r —>■ cy) will always be close to this number and consequently 
at least two orders of magnitude below the future limits (between 10~^ and 10~^) which 
are currently foreseen. 

To be more precise, the branching ratios as defined in Eq. (fTOl) have to be evaluated in the 
basis in which the charged leptons and the heavy Majorana neutrinos are diagonal. In this 
case Md has to be replaced with Md = Md Ui, where Mr = Mf"^ V^. Obviously, 
Vr drops out of M\j Md and does not influence LFV. Now consider again the case that 
Ui is non-trivial and given by a 23-rotation, which we showed above to keep the scaling 
predictions of = Ue3 = and to change c to c given in Eq. ([6]). One easily finds that 
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iMj,Mn)2i 



cos + c sin 9. 



23 



c cos6' 



23 



sin 9. 



23 



(12) 



This expression is nothing but 1/c^ and therefore we recover the relation Eq. (fTTj) between 
the "double ratio" BR(r e7)/BR(/i cj) and atmospheric neutrino mixing. 



^Note that there is no need for Mn to obey scahng because it would be singular in this case, invalidating 
simple seesaw. 
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It is interesting to ask whether scahng can be apphed to non-standard neutrino scenarios. In 
particular, the possibility of neutrinos being Dirac neutrinos and the presence of additional 
light sterile neutrino species is discussed frequently. 

Dirac neutrinos could be accommodated by scaling if the neutrino mass matrix would take 
the form given in Eq. ([9]). The resulting neutrino oscillation phenomenology (obtained 
by diagonalizing MjyMi)) would be again an inverted hierarchy with tan^^23 = and 
Ue3 = = 0. The reason is simply because the resulting mass matrix possesses an 
eigenvalue with a corresponding eigenvector of (0, — 1/c, 1)-^. 

Sterile neutrinos would require enlarging Ai^, from being a 3x3 matrix to a {3+ns) x {S+Us) 
matrix, where is the number of additional sterile neutrinos. The recent results of the 
MiniBooNE experiment [H] seem to indicate that > 2 [15j. We can modify the scaling 
condition from Eq. ([1]) to include also the sterile neutrinos: 

^e/i ^/X/i f^Tfl ^SlfJ, '^S2lJ. ^ (13) 

TTIq-j- T^^jlt TfLq-j- TfLg-^j- Tn^^T 

The result is a mass matrix with a zero eigenvalue having an eigenvector (0, —1/c, 1, 0, 0, . . .)'^ 
The mixing scenario is described by f/es = 0, lU^s/Ursl"^ = 1/c^ and Ug^^s = Ug^s = . . . = 0. 
Leptonic CP violation is important in order to allow scenarios with two sterile neutrinos 
to survive constraints from current data [15j. In total there are five CP phases in this 
case, and only three of them are unphysical due to the vanishing mixing matrix elements. 
For two sterile neutrinos, the resulting scenario would correspond to an inverted hierarchy 
of the three mostly active neutrinos and two heavier, mostly sterile neutrinos. It is in 
the language of Ref. [TB] the scenario "SSI" , which has from all possible mass orderings 
the smallest predictions for the various mass-related observables (neutrinoless double beta 
decay, neutrino mass in KATRIN and the sum of masses in cosmology) [16]. 



3 Scaling Predictions for Astrophysical Neutrinos 

It has recently been recognized that measuring flux ratios of high energy astrophysical 
neutrinos [17] is an alternative method to determine neutrino mixing parameters fiSi [T9l 
[20] . In particular, one expects from astrophysical pp or p~f processes, which generate 
pions and kaons, an initial flux composition of the form $g : $° : $° = 1 : 2 : 0, where 
$1^ with a = e,fi,T is the flux of neutrinos with flavor a. Neutrino mixing modifies 
the flavor composition and in terrestrial neutrino telescopes such as IceCube [21] one can 
then measure flux ratios and thereby obtain information on the neutrino parameters. The 
measurable neutrino flux is given by 

$a = E^"/5'^° withp,^ = 5]if/„,n[/^,|2 . (14) 

In the limit of maximal atmospheric neutrino mixing and vanishing Ues, the composition 
1 : 2 : is transformed into 1:1:1, independent of the solar neutrino mixing angle. Small 
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deviations from 6*23 = vr/4 and 613 = lead to [HI [20 



<l>e:$^:$r = l + 2A:l-A:l-A, 

1 1 /I \ (15) 

where A = - sin4^^i2 jt/esl cos 5 H — sin^ 2^^i2 ( sin^ ^23 ) ■ 

4 2 V 2 / 

Thus, there is a universal first order correction in terms of the small parameters |?7e3| and 
I — sin^ ^23- With the current la {3a) ranges of the oscillation parameters one finds that 
-0.036 < A < 0.057 (-0.097 < A < 0.113). The ratio of electron neutrinos to the other 
flavors is therefore a probe of ^23, ^13 and cos S. Note that in the definition of A in Eq. ( |T5|) 
the factor in front of | — sin^ ^23 is larger and has a smaller range than the one in front 
of |?7e3| cos 5. To be precise, for the allowed 3a range of solar neutrino mixing, ^ sin 4^12 
ranges from 0.12 to 0.21, whereas |sin^26'i2 ranges from 0.38 to 0.48. Consequently [20] . 
the sensitivity to deviations from maximal atmospheric neutrino mixing is better than the 
sensitivity to deviations from \Ue3\ = 0, which in addition gets smeared by the dependence 
on the CP phase 6. 

We conclude that scaling - predicting only non-maximal 6'23 and neither 6 nor ^^13 - will 
have particularly simple, interesting and potentially testable phenomenology at neutrino 
telescopes. To go into more detail, let us introduce the small parameter 

6 = ^-^23, (16) 

which can be linked to the scaling parameter c via l(c^ — l)/(c^ + 1) = 1 — sin^ ^^23 — £• 
The result for the flux ratios is 

$e : : $r = 

1 + 2Ci2 S12 (^23 ~ '^23) • 2 (1 — 2C23 '^23 ~ '^12 ^12 (^23 ~ '^23) ^23) • ^23 (l + (1 ~ '^12 '5l2)(c23 ~ ^23)) 

~ 1 + 4 c?2 e : 1 - 2 2 e : 1 - 2 cj^ s?2 e , 

(17) 

where we have given the exact expression and the expansion in terms of e to first order. 
The electron neutrino flux $e receives no quadratic correction and the second order term 
for $„ is 4 (1 — c^2 ^12) while that for is identical but with opposite sign. In Fig. [T] 
we plot - using the exact probabilities - the flux ratios ^e/^fi and $^/$tot5 where $tot 
is the total neutrino flux, as a function of the scaling parameter c. For c = 1 we obtain 
^e/^^i = 1 and $^/$tot = |- Quite large deviations from these values are allowed, and the 
dependence on solar neutrino mixing is very weak. 

Finally, we note that if there are sufficiently fast non-standard decay modes of the neutrinos, 
and only the lightest state Ui {i = 1 for normal ordering, i = 3 for inverted ordering) survives 
and is detected, the fluxes obey the relation [22] '■ '■ = \Uei\^ '■ |f^/tiP : l^riP- In 



- • ^ T I ^ 62 I • \^ fit] 

case of scaling (inverted hierarchy and ^13 = 0), this simplifies to 

<l>e : $^ : $r = : sin^ ^23 : cos^ ^23 = : 1 : . (18) 

There are no electron neutrinos and the ratio of muon to tau neutrinos is tan^ ^23 = l/c^. 
We plot for the case of decaying neutrinos in Fig. [1] the ratio of muon neutrinos to the 
total fiux, which is simply equal to sin^ ^23- 
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Figure 1: Flux ratios of high energy astrophysical neutrinos as a function of the scahng 
parameter c. The red (sohd) hues are ^e/^fi for the standard case of an initial 1:2:0 
composition, the green (dashed) hues are $tot in the standard case and the violet (dot- 
dashed) line is for neutrino decay. The first two cases are plotted for the best-fit values 
and the upper and lower 3cr values of the solar neutrino mixing angle 6i2- 

4 Non-maximal atmospheric Neutrino Mixing and van- 
ishing Ue3 from other Scenarios 

The question arises if we can obtain the values 6*13 = and 623 7^ 7r/4 by breaking or 
modifying scenarios with initial ^13 = and 623 = n/A. Explicit breaking of the symmetry 
in a mass matrix, RG effects and contributions from the charged lepton sector are appealing 
possibilities, which we will now comment on. 

Considering first radiative corrections, we note here that 613 and 623 are inversely propor- 
tional to Am\ [8] , and due to this one would expect that in general RG corrections to them 
are of the same order and in addition small. Therefore, if initially ^13 = ^^23 — 7r/4 = 0, then 
keeping at low energy ^13 zero but having simultaneously 623 non-maximal would require 
rather special values of the other parameters and unnatural cancellations in particular 
between the CP phases. Hence, 623 7^ vr/4 and 613 = should hold initially and we end 
up again with the question of how these pecuhar values arose, which leads us back to the 
scaling hypothesis. 

Turning to explicit breaking of the symmetry in a mass matrix leading to ^13 = and 
^23 = 7r/4 requires a glance at /i-r symmetry ^23|, which imprints the following form on 
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the mass matrix: 



/ A 



B 



M 



u = mo 



B 



D 



E 



(19) 



E 



Maximal 6*23 and zero Ues are predicted, but only if the eigenvalue belonging to the eigen- 
vector (0, —1, 1)"^ is the largest or smallest one, which would then correspond to the normal 
or inverted mass ordering. In case of scaling, the eigenvector (0, — 1/c, 1)"^ belongs automat- 
ically to the zero eigenvalue, therefore there is no such ambiguity. For a normal hierarchy 
of the neutrino masses the parameters in Eq. (fT9l) have to fulfill D, E ^ A, B. If in this 
case the fi-r symmetry is broken such that only the fifi and fir entries differ, but the e/i 
and er elements stay identical, then it turns out that Ues is small but non-zero (to be pre- 



cise, it is of order AttIq/ Am\), while 6*23 deviates from maximal by roughly AttIq/ Am\ 

[21] . This way of breaking has no analogue for an inverted hierarchy or quasi-degenerate 
neutrinos, for which rather tuned breaking scenarios are required in order to end up with 
^13 = but ^23 ^ vr/4. 

Contributions from the charged leptons arise when A^j, is /i-r symmetric and in the limit 
of a diagonal charged lepton mass matrix 6*13 = and ^23 = 7r/4 would result from U = 
Ul = Uy. It is commonly assumed that Ui contains only small angles. Introducing for 
the sines of these angles the abbreviations sin6'fj = Ajj, one finds in first order of these 
small parameters that [25|, [26] 



where 0i and 02 are CP phases. Consequently, 6^13 = and 6*23 7^ vr/4 would require del- 
icate interplay of the angles and phases in and would lead to a rather unnatural form 
for it. In particular, for the natural case of Ui being CKM-like, i.e., A12 S> Ai3^23, the result 
6*13 = and 6*23 7^ vr/4 can not be achieved. Similar statements hold for the opposite case 
in which in the limit Uy = 1 the charged lepton sector would suffice to generate 61^ = 
and 023 = 7r/4 mU = U} = U} 

We conclude that it seems rather unnatural to obtain [/es = and 6'23 7^ 7r/4 within broken 
/x-r symmetry. These predictions for the mixing angles also occur in models based on the 
flavor symmetry Le — L^^ — Lr- However, as we will show in the next Section, various tuning 
problems show up for this Ansatz. 

5 Inverted Hierarchy: Scaling vs. the L^ — L^ — Fla- 
vor Symmetry 

Note that scaling is an Ansatz for the inverted hierarchy which is fundamentally different 
from the flavor symmetry — — Lr [6], which is usually "blamed" for it. A detailed 




|t/e3|-75 Ai2-Ai3e*^i , 
sin^ 6*23 ^ I + A23 cos 02 - I (A12 - A13) + \ cos 0i A12 Ai3 , 



(20) 
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comparison is therefore a worthy exercise. A mass matrix obeying the flavor symmetry 



Mu = mo 



I COS 6^23 -sin 6123 \ 

cos^ 
\ - sin 6*23 y 



where mo = JAmj^ , (21) 



and generates = Ue3 = as well as non-maximal atmospheric neutrino mixing given 
by ^23- However, it also predicts maximal solar neutrino mixing and vanishing Amg = 
m2 — mf, which is in contradiction to observation. Therefore, in contrast to scaling, the 
symmetry needs to be broken to achieve correct phenomenology, which imposes three 
problems [21 



(i) the breaking terms in the mass matrix have to have at least 30% the magnitude of 
the terms allowed by the symmetry. The reason is that in the inverted hierarchy the 
ee entry of A4,y (the effective mass) is required to be larger thar|§ J Am\ cos 2^12; 



(ii) the (large) breaking of the symmetry in the mass matrix is always connected with 
fine-tuning because usually the required large deviation from maximal du is con- 
nected with the small ratio of the solar and atmospheric mass-squared differences. 
For instance, let us add to the matrix in Eq. fl?I]) the following (/i-r symmetric) 
perturbation: 

f a b b 
e b d e 
\ b e d 

If for simplicity we set 623 = 7r/4 in Eq. fl?T]) then maximal atmospheric neutrino 
mixing and Ue3 = will remain, but the ratio of mass-squared differences is now 
AmQ/Am\ ~ ^/2(a + d + e) e while solar neutrino mixing is governed by sin^^i2 — 
l/^/2 -{a-d-e) e/8. Hence, if one wants to reproduce the best-fit values of the 
oscillation parameters in Eq. ([5]), the fine-tuned condition {a + d + e) / {a — d — e) ~ 
0.0027 has to be fulfilled; 

(iii) if the symmetry is broken by contributions from the charged lepton sector (note 
that in this case in addition breaking in the neutrino sector is necessary to generate 
non- vanishing AttIq) a CP violating phase appears in the expression for the now 
non-maximal solar neutrino mixing angle, which is required to be close to zero: using 
again the natural choice of a CKM-like Ui (see also Eq. fl20|) ) leads to the formula 



sin^ 6^12 ^ ^ - cos cos 6^23 A12 , (22) 



^The fact that L^. — — predicts an effective mass near the lower end of its allowed range, whereas 
scaling admits (m) to take any of its allowed values, represents a possibility to distinguish the two Ansatze. 
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where A12 is the leading 12-rotation in Ui. From the experimentally observed sin^ 612 ~ 
0.3 it follows for the natural value of A12 — 0.2 that cos0 has to be tuned to be very 
close to one. Leptonic CP violation in oscillation experiments is proportional to sin (p 
[271 [26] and very much suppressed. 

All these fine-tuning problems occurring in Le — — Lr do not occur in scaling, which 
therefore represents a presumably better Ansatz for the inverted hierarchy. Moreover, as 
we will elaborate upon in the next Section, scaling can easily be obtained in models based 
on discrete flavor symmetries, which are currently intensively studied [28]. 



6 A Model for Scaling and Phenomenological Conse- 
quences 

We consider now a seesaw model based on the D4 x Z2 flavor symmetry which was proposed 
in Ref. [2] to generate scaling. The particle content together with the quantum numbers 
under D4 x Z2 is shown in Table [H The superscripts +, — refer to the transformation 
under Z2 and the rest are the 1)4 representations. For the mathematical details of the D4 
group see for instance Ref. [29]. Apart from the usual Majorana neutrinos N^^^^^r, the right- 
handed charged leptons cr, ^ir, tr and the lepton doublets Le^^^r, one has to introduce five 
Higgs doublets 0i, 2,3,4,5- In the Appendix we show as a proof of principle that the D4 x Z2- 
invariant Higgs potential can be minimized with Higgs masses having values above current 
limits. From the assignment in Tabled] the following Lagrangian is obtained: 

^ = kieR^4>i Le + JlR{k24>i - hh) Lf, + T]^{k2(f)i + ^(1)^)1^ 

+ iVe 01 Le + /l2 A^02 + /is (04 + 05 ^/.) (23) 

+i {NlK Mi + N;N^^M2+ KN^ Ms) + h.c. 
Hence, the neutrino Dirac mass matrix can be written as 

/ ae^^ \ 

Md= b d e t;„k , (24) 

y y 

and both the charged lepton and Majorana mass matrix are diagonajf]. Note that multi- 
Higgs models as the one analyzed here typically predict flavor changing neutral currents 
and LFV in the charged lepton sector at dangerous levels. Here the model has a diagonal 
charged lepton mass matrix which renders processes like fi —>■ ej suppressed either by the 
GIM mechanism or by the masses of the heavy right-handed neutrinos. Note that the 
model presented here is non-supersymmetric. A supersymmetric version would have LFV 



■^There is another model based on D4 x Z2 presented in Ref. [2 . The heavy Majorana mass matrix 
is arbitrary in this case, which has therefore httle predictivity for leptogenesis. 
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Field 


D4 X Z2 quantum number 


Le 


1^ 


Cr, Ne, 01 


ir 


iV^, 02 


1+ 

-■-2 




I2" 




I4" 




2+ 








2- 







Table 1: Transformation properties under D4 x Z2 of the particle content of the model. 



via off-diagonal slepton mass matrices generated by the mechanism described in Section 
[2] (see Eq. (ITUl) ). Because from Eq. is a special case (recall that e = d/c) of the 
form given in Eq. (Q one would in this case obtain the characteristic relation between the 
LEV charged lepton decays and atmospheric neutrino mixing from Eq. ffTTl) . 
In Md shown in Eq. flM|) we have already included one complex phase. It is easy to show 
that with diagonal charged lepton and Majorana mass matrices there is only one complex 
phase in the model. Using all this, we can calculate the neutrino mass matrix using the 
type I seesaw formula to obtain 



-w. 



( M2 

Ml 



V 



hd 
he 



2 hd he 
d^ de 
de e^ 



\ 



I 



(25) 



Note that the third heavy neutrino mass M3 does not appear in M.v, i-e., effectively we 
are dealing with a 2 x 3 seesaw. The low energy mass matrix apparently obeys scaling 
with c = d/e. We therefore have = Ues = and tan^ ^23 = e^/d^. By evaluating 
the rephasing invariant Im {meerrirTrnl^m*^} [30j both with Eq. fl25|) and with the usual 
parametrization of the PMNS matrix, one finds a compact relation between the parameters 
in Md and Mr and the low energy observables: 



1 A 2 
- Ami 

4 ^ 



A 



.m|, sin^ 023 sin^ 26*12 sin 2a = 




aH^ e' 



sin2(y9 



(26) 



The effective mass governing neutrinoless double beta decay is 



(m) 



Mo 



M 



Ml 



'-a^e^'^ + h'' 



Am A V 1 - sin^ 26*12 sin^ a 



(27) 



Let us focus now on high energy phenomenology in terms of leptogenesis [3T] . The decay 
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asymmetries of the heavy neutrinos A^i,2,3 into final states with fiavor a = e, fi,T are [32] 

T{N,^<j)Q-T{N,-^<f)U^) 



/ 1 fl + X 

where f(x) = \fx 1 H (1 + x) In 

V 1 — X \ X 

With the very restricted form of the Dirac mass matrix from Eq. fl2^ it turns out that 
only two of the nine possible are non-zero. Those are 



e\ = ^\? sin 2^ i {Mil Ml) ^ sin2v. ^ 



^ .2^W2^ 2 sin2¥P f{Ml Ml) ~ — , sm2¥; -i 



(29) 



where we also gave the limits for M| ^ Mf. Again the third heavy neutrino does not play 
any role. The leptogenesis phase ip is identical to the low energy Majorana phase a, so that 
we can expect a correlation between the baryon asymmetry and the effective mass governing 
neutrinoless double beta decay. Let us focus in the following on the very typical case that 
El = el governs the baryon asymmetry. We need to specify the magnitude of the heavy 
neutrino masses. The common mass scale v^y./M2 in the light neutrino mass matrix from 

Eq. (|25ll should be Am\ ~ 0.05 eV, which brings us to the choice M2 = (5- 10^^ ^5- 10^^) 
GeV, when the entries of the mass matrix lie between 0.1 and 10. We further note that in 
case of an inverted hierarchy all entries of Ai^ are of the same order of magnitude, so the 
ratio M2/M1 should not be too large. Nevertheless, there should be a moderate hierarchy 
in order to avoid the complications and tuning issues of heavy neutrinos close in mass, 
so for definiteness we choose Mi = (| ^ |) M2. This in turn means that fiavor effects in 
leptogenesis [32j do not play a role. In this limit, we can estimate the baryon asymmetry 
as 

Yb ^ csp ^'^("^i) ' (30) 

g 

where g* = 122.75, csp = —44/87 and k can be parameterized as [33] 

, , /0.55-10-3eVV''' 

f^{mi) = z 

\ mi J 

with rhi = (M^) M|))ii/Mi. In Fig. [2] we show the results of an analysis in which we 
searched with the ranges of Mi and M2 specified above for values of a, b, d, e and which 
generate the correct neutrino mixing phenomenology as defined in Eq. ([5]). Having found 
such parameters we evaluate the baryon asymmetry, which should lie in the range (8-^10) ■ 
10~^^ [31]. We see that there is as expected a correlation between Yb and (m) and that 
several of the points generate the correct baryon asymmetry, both in sign and magnitude. 
The particular choice of the parameters demands for successful leptogenesis the effective 
mass to lie around its maximal allowed value, (m) ^ 0.045 eV. 
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Figure 2: Scatter plot of the effective mass against the baryon asymmetry of the Universe 
in the model from Section El 

7 Summary 

In summary, we present a detailed investigation of the hypothesis that the Majorana neu- 
trino mass matrix obeys a scaling law as a way to understand current neutrino observations. 
Two consequences of this hypothesis are that (i) the neutrino mass ordering is inverted 
rather than normal and (ii) both U^s, and the lightest neutrino mass vanish. These results 
are invariant under renormalization group extrapolation and are therefore stable under ra- 
diative corrections, which distinguishes the scaling proposal from many others motivated 
by family symmetries or texture zeros. The effective mass governing neutrinoless double 



beta decay can be as large as \J Am\ ^ 0.06 eV. Another distinguishing prediction of 

scaling is that the value of the atmospheric mixing angle is not necessarily maximal even 
though Ue3 = 0. This is in contrast to models with /x-r symmetry, which provide a simple 
way to understand both maximal atmospheric mixing with a very small Ues- Models with 
approximate or broken /i-r symmetry always correlate non-vanishing of Ues with deviations 
from maximal 623. We note that recent analyzes of the available neutrino data tend to 
favor non-maximal atmospheric neutrino mixing. We also compare the scaling hypothesis 
to the L^ — L^ — Lj. flavor symmetry which is very often utilized to understand an inverted 
hierarchy. Unlike the scaling Ansatz, fitting observations requires a fine-tuning of mass 
matrix elements/perturbation parameters. An interesting aspect of our hypothesis that it 
is invariant under any possible rotation of the basis in /x-r space, e.g., coming from the 
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charged lepton sector. This may make it easier to construct models that obey scahng. We 
also note ways to test scaling using high energy astrophysical neutrino fluxes. We discuss 
a particular class of seesaw models based on the D4 x Z2 flavor symmetry that realize 
the scaling hypothesis. There are (i) no dangerous flavor changing neutral currents in the 
lepton sector, (ii) the leptogenesis phase is identical to the low energy Majorana phase and 
(iii) successful leptogenesis requires the effective mass in neutrinoless double beta decay to 
be larger than 45 meV. 
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A Scalar Potential 

The model we study in Section E] is defined by the transformation properties given in 
Table [H It includes five Higgs doublets 0i, 02, 03, 04 and 05. We will show now that 
the potential involving all these fields can be minimized with a realistic vev configuration 
leading to a phenomenologically viable spectrum of scalar bosons below mass bounds from 
direct production. Five Higgs fields correspond to ten charged and ten neutral ones, or to 
20 physical degrees of freedom. Of the ten charged degrees of freedom, two are eaten by the 
electroweak W bosons, leaving eight degrees of freedom, i.e., four complex charged scalars. 
Of the ten neutral degrees of freedom five are even under charge conjugation. These give 
five real physical scalars. The other five neutral degrees of freedom are odd under charge 
conjugation. One of these is eaten by the electroweak Z boson, leaving four physical 
pseudoscalars. With the assignment given in Table [T] the most general x Z2-invariant 
potential with real coefficients reads: 

3 3 
1=1 1=1 

+ Ai2(0l0l)(0502) + Ai3(0l0l)(0j03) + A23 (0^02) (0^03) 
3 

+ ^Ki(0j0i)(0l04 + 0^05) +ai[(0l02)^ + /i.e.] +a2|0l02r (Al) 
i=l 

+ asKcplcpsf + h.c] + 0410203^ + "5(0405 + 0504)^ + a6(0405 " 0504)^ 

+ 07(0404 - 0^505)^ + a8(0404 " 0^505) (0103 + h.C.) + a9[(0504)^ + (0^05)^ + h.c] 

+ aio(|0204|' + 10^051') 
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+ «ll[(0l04)' + {4<P5? + h.c] + ai2(|0l0, " 

+ ai3[(0304)' + {4<p5f + h.c] + au{\(f)l<pAr + 



41 + 
2 



4051') 



This potential can for instance be minimized by choosing the VEV configuration: 



= y j for i = 1,2,4,5 and (^a) = 
if one uses the following numerical values for the coefficients: 



Ai = 2.62904 
A4 = 2.93754 
A23 = -0.455254 
K3 = -0.778689 
03 = -2.649 
ae = -1.82076 
ag = -0.853325 
= 0.68049 
ai5 = 1.70881 



2 

/^2 



0.661705 w 



wk 



A2 = 2.91805 

^2 

Ki = 



2 — ^.^^^.^^ , 

Ai2 = 2.94576; 
0.190918; 
-1.58251; 
04 = 2.32869 ; 
aj = -0.097017; 
aio = 0.393294; 
ai3 = -1.29704; 
ai6 = -0.322439 ; 
-0.971516 



wk ' 



A3 = 
Al3 
K2 = 

a2 = 

"5 = 
"8 = 

an 

"14 

fil = 
2 

/i4 = 



2.53936; 
= -1.29472; 

2.25757; 

-0.77218; 

0.690183; 

-1.38637; 
= -0.404394 ; 
= 0.764254 ; 



(A2) 



(A3) 



0.619434 1;, 



wk ' 



1.4775 V 



wk 



One finds masses of 552 GeV, 390 GeV, 362 GeV, 240 GeV and 147 GeV for the scalars, 
682 GeV, 607 GeV, 460 GeV and 322 GeV for the pseudoscalars and 404 GeV, 363 GeV, 
273 GeV and 155 GeV for the charged scalars. These scalars are in general superpositions 
of all five Higgs fields, except for the 240 GeV scalar, which is only a linear combination 
of 04 and 05. Note that our model is only focussing on the lepton sector. Confronting 
the obtained Higgs mass values with limits stemming from rare meson decays would mean 
to construct a full model including also the quark sector and carefully perform a lengthy 
study of the diagrams leading to fiavor changing neutral currents taking into account all 
five Higgs doublets. Within our model the specific choice of Eq. (lA2p leads in the Dirac 
mass matrix from Eq. (p^ to c? = e and therefore maximal atmospheric neutrino mixing. 
Other values are of course possible, which would lead to slightly different scalar masses 
and parameters in Eq. (lA3p . 
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